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Abstract 

We work in the framework of a simple vector-tensor theory. The parametrized post-Newtonian 
approximation of this theory is identical to that of general relativity. Our attention is focused on 
cosmology. In an homogeneous isotropic universe, it is proved that the energy density, /9^, of the 
vector field A, and its pressure, p^, do not depend on time, and also that the equation of state is 
= —Pjy- This means that, in the theory under consideration, there is a cosmological constant, 
which is not vacuum energy, but the dark energy of the cosmic vector field A, whose evolution is 
classical. 
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I. INTRODUCTION 



We begin with our motivations to study vector-tensor theories. As it is well known, the 
analysis of the WMAP (Wilkinson microwave anisotropy probe) data has pointed out a set 
of anomalies in the cosmic microwave background (CMB) angular power spectrum, e.g., the 
planar character of the octopole and its alignment with the quadrupole, 1^2], [s], 0, [5] the 
asymmetry between the north and south ecliptic hemispheres jo, Q, Is, 9, lo[, and so on. 
These anomalies strongly suggest that the distribution of CMB temperatures deviates from 
statistical isotropy at very large angular scales. Most of these anomalies have been recently 
explained [111 by using divergenceless peculiar velocities (superimpositions of the so-called 
vector modes) in the framework of general relativity (GR); nevertheless, in this theory, 
the angular velocity of these divergenceless motions decays during the matter dominated 



era and, consequently, such motions should be generated a 



ose to z = 1100 to produce 



CMB anisotropy after decoupling and before decaying [11|, Il2|. This condition seems to be 
necessary to make possible the explanation of the low i CMB anomalies proposed in 11| ; 
nevertheless, no specific processes are known to create the required divergenceless velocity 
fields around z = 1100. By this reason, we propose the use of some suitable theory involving 
more vector modes than GR. In this alternative theory, the evolution of the angular velocity 
could be very different and more appropriate to explain anomalies. These comments suggest 
the choice of some vector-tensor theory, whose four-vector field leads to new vector modes 
associated to its spatial part. According to these ideas, we are performing a study of the 
vector-tensor theories obtained from the Lagrangian described in the next section 13|. Other 
theories will be considered in the future. Our first outcomes on the evolution of vector modes 
in vector-tensor theories were presented in fl^ . 

Before studying the evolution of the vector modes in a given vector-tensor theory, 
we should explain some well known observations; e.g., it should have an appropriate 
parametrized post-Newtonian (PPN) limit compatible with solar system observations. In a 
first step, we are studying some vector-tensor theories leading to the same PPN parameters 
as GR. Moreover, in the chosen theory, the homogeneous and isotropic cosmological back- 
ground should explain, at least, all the observations already explained by the concordance 
model. The main goal of this paper is the description of a theory with these characteristics. 
We have found one of them, in which the vector field plays the role of the cosmological con- 



2 



stant. No vacuum energy is then needed to explain SN la observations, CMB anisotropies, 
baryonic acoustic oscillations, and so on. These surprising results have motivated this brief 
paper. 

We think that the new classical origin of the cosmological constant is a very interesting 
result. It deserves attention even if CMB anomalies cannot ever be explained in the frame- 
work of the selected vector-tensor theory. The evolution of the vector modes (perturbations 
of the background four-vector used here) is not required in this paper. It is necessary to 
try to explain anomalies (our initial motivation), but this possible explanation is yet to be 
studied. Results will be presented elsewhere. 

In this paper, G, a, r, and z stand for the gravitation constant, the scale factor, the 
conformal time, and the redshift, respectively. Greek (Latin) indexes run from to 3 (1 
to 3). Whatever function ^ may be, ^* is its present value and ^ {^') denotes its partial 
derivatives with respect to r (z). Quantity pc is the critical density, pr (pm) is the radiation 
(matter) energy density, and pb = Pr + Pm is the total background energy density of the 
cosmological fluid. Since we work in a fiat universe, the present value of the scale factor a* 
is assumed to be unity and, then a = (1 + z)'-^. Quantity w is the ratio between pressure 
p and density p. Units are defined in such a way that the speed of light is c = 1. Finally, 
indexes are raised and lowered with the space-time metric. 



II. VECTOR-TENSOR THEORIES 



In GR, there is a tensor field f?^^ which plays the role of the space-time metric. In the 
vector-tensor theories, there are two fields, the metric q,.^ and a four- vector A^. Several of 

n n 

these theories have been proposed (seejl3|, and references cited there). We have first 



considered those based on the action 
I = (IGn 



13|: 



Gy' J{R + CoA^A^R + r/i?^,AM^ - eF^.F^'^ + 7 V.A^VM'^ + L^) ^ d' 



X 



where u, fj, e, and 7 are arbitrary parameters, R, R^u-, 9, and are the scalar curvature, 
the Ricci tensor, the determinant of the g^i, matrix, and the matter Lagrangian, respectively. 
The symbol V stands for the covariant derivative, and F^^, = V^A^, — Vi^A^. In action ([T]), 
it is implicitly assumed that the coupling between the matter fields and is negligible. 



A more general action is given in 



Iq . It involves a new term of the form X{A^A^ + 1), 



where A is a Lagrange multiplier. With the help of this term, the vector A^^ is constrained 
to be timelike with the unit norm. From this last action, the field equations of the so- 
called Einstein-Aether constrained theories can be easily obtained. The applications of 



these theories to cosmology are discussed, e.g., in 
is used in 



161 ]. A mass term of the form irP'^A^A^ 



m 



17| to explain cosmic acceleration with a massive vector field. The same is done 



for an unconstrained theory based on action ([T]). Recentl y, o ther theories involving 



vector fields have also been applied to cosmology (see, e.g., |l9l . l20|]). Fortunately, among 
all the proposed vector-tensor theories, we have found a simple unconstrained one, based on 
the action ([1]), which leads to a new interpretation of the concordance model in the absence 
of vacuum energy. 

It can be easily proved that, for = 0, = 7, and arbitrary e, the PPN parameters of 



the theory based on the action ([T]) are identical to those of GR 



13| 



the free parameters we easily get the following field equations (see 13|): 



2l|. For this choice of 



1 



G 



R 



pi/ 



SttGT, 



(2) 
(3) 



where Gfj,u is the Einstein tensor^ 



Q^t^u , Q^i?u , and Gjfj is given in [13|]. Taking into account the identity 



eO^fj, and the explicit form of 



Eq. ([2]) is easily rewritten in the form: 



fjV^iV ■A) + {2e-fi)V'F^, = 



(4) 



(5) 



where V ■ A = VqA". By using the identity (jlj) and Eq. ([5]), the energy-momentum tensor 



of the field A^; namely, the tensor = T^/SttG can be written as follows: 



- |(V ■ + {2e - rj) (a^V^F^/j - -F^/^F"^ 
i2e - r])g^''{F^aF,p + A^V pF,^ + A,V pF^^) 



(6) 



where 77 = rj/SnG, and e = e/SnG. Finally, in the cosmological the energy- 

momentum tensor of the cosmological fiuid, which involves both matter and radiation. It is 
worthwhile to notice that the relation 



VT^ = 



(7) 



is satisfied fsee 



131 ]) and, consequently, taking into account Eq. ([3]) and the identity V^G 



0, the relation V^T^jy = is also satisfied. This means that matter and radiation evolve as 
in the standard Friedmann-Robertson- Walker model of GR and, as it is well known, after 



22l|). the following equations are 



the annihilation, which took place at 2; ~ 10^'' (see 
valid: 3pr = pr = p*(l + z)^, pm = Pm(l + zf, and Pm ^ 0. 

For u) = f/ = 7 = and e ^ 0, the term —sF^yF^^^ involved in the action ([1]) has the 
same form as the term appearing in the action of an electromagnetic field in the absence of 
sources. In this case, Eq. ([5]) reduces to 'S/^F^^ = and the energy- momentum tensor of the 
field Af, is T^t = 2£[^°^F^^F^/3 - (1/4)^^^F„^F"^]. As it is well known, the field is called 
the potential vector of F^^. The equations of the theory can be written in terms of the field 
Ff^i, and its derivatives and, consequently, classical theoretical predictions and observations 
refer only to the field F^^^. Some indetermination in the potential field is unavoidable (the 
potential vector can be arbitrarily fixed by using appropriate gauge conditions). Similarly, 
for other values of the parameters involved in the action ([T]), the vector field may also 
play the role of a potential generating other measurable physical fields involved in the field 
equations. The role of vector and the indetermination in it depend on the properties of 
the Lagrangian (invariance under transformations of A^). For example, in the case a) = 0, 
r/ = 7, and 2e = fj 0, the energy-momentum tensor reduces to = — ^(V ■ A)^g^i, 
and Eq. ([5]) gives V^(V - A) = 0. According to this last equation V ■ A is constant and, 
consequently, tensor T^^, has the same form as the energy-momentum tensor corresponding 
to vacuum; namely, one has T^^ = —p^g^u, where p^ = f (V ■ A)"^ = constant 7^ 0. This 
means that the resulting theory is equivalent to GR plus a cosmological constant. In this 
theory, is not a classical field to be determined either with the field equations or with 
observations, it is a field playing the role of a potential vector for the scalar field V ■ A, which 
is the only field (apart from the metric) involved into the field equations of the theory. 

In this paper, we pay particular attention on the theory corresponding to u) = 0, = 7 7^ 
0, and ^ 2s ^ fj. For this choice of the free parameters of the action ([T]), the field equations 
(see above) are more complicated than in the simple case 2e = fj; for example, according 
to Eq. ([5]), the scalar V ■ A is not constant, excepting some special physical systems as, 
e.g., a homogeneous and isotropic background universe. The evolution of the cosmological 
background is studied in the next sections. 
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III. COSMOLOGY: BASIC EQUATIONS AND COSMOLOGICAL CONSTANT 



In this section, we focus our attention on a homogeneous and isotropic cosmological 
background. In the flat case, the line element is 

ds^ = a?{-dT'^ + Sijdx'dx^) . (8) 

Moreover, homogeneity and isotropy require a vector field whose covariant components are 
(Ao(r), 0, 0, 0). Then, tensor F^^, vanishes and 

V-A = --(a, + 2-A^ ■ (9) 



a? \ a 



hence, Eq. ^ reduces to the relation (i(V ■ A)/dT = 0, which can be rewritten in the form: 



• 2 

a a 



a a? 



Ao + 2Ao - - 3- = . (10) 



Equation ([6]) allows us to find the components of T^. By using these components, Eqs. ([3]) 
can be written as follows: 



and 



where 



3— = 87rGa2(p, + pJ (n) 



-2- + ^ = 87rG'a2(p,+pJ , (12) 



2 



P. = -p. = |(V.^f = ^(A„ + 2^A„j ; (13) 

hence, = —1, as it occurs in the case of the cosmological constant {wa = —1). 

Since we have found that — in the cosmological case — the energy-momentum tensor of 
the field A has the form 

T;t = -pAQp^y ' (14) 

Eq. ([7]) implies the relation = 0; hence, the energy density is constant. This is in 
agreement with the relations p^ = ^(V ■ A)"^ and d(V ■ A)/dT = previously obtained. 

Evidently, if the energy density of the evolving field A is constant and the relation p^ = 
—p^ is satisfied, we can state that this field acts as a cosmological constant. The possible 
values of p^ are discussed in the next section. 
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IV. SOLVING THE BASIC COSMOLOGICAL EQUATIONS 



In order to integrate Eqs. (fT0|) - (fT2|) . the following new variables are appropriate: ( = 1+z, 
Hi = t' , 1/2 = Aq, and 1/3 = A'q. In terms of these variables, the energy density is 



2 



and Eq. ( |TT1) can be rewritten as follows: 



Pa = ^(i/3C-22/2)', (15) 



vC\y,C - ly^f + 2yie {CpI + P™) = ^ , (16) 



whereas Eq. f[T^ leads to 



y[ = -|i {47rGC^ [2ylpX' - Sr^iy^C - 2^2)'] + 9} . (17) 



Finally, the second order differential equation (fTOj) leads the following system of first order 
equations: 

1/2 = ys (18) 



and 



/ _ 2?/2i/i + y'lC {ysC - 2?/2 
yiC 

From Eqs. f|T6l) and f|T71) one easily finds 



/ _ ^y2yi -r j/jc, yyai-, - -^^2; /.qx 
^3 — :7T5 v-'-^J 



?/? 3 



which can be easily integrated to get 



yi 



(4pX' + 3p;,C') , (20) 



1/2 

(21) 



.87rG(pX' + p^e + C) 

By combining Eqs. (fT5|) . (fT6l) . and (I2T1) we can easily get the equality C = p^; therefore, 
taking into account that C is an arbitrary integration constant, the constant can take on 
any value and, in particular, we can fix the value Pa = P*^^ = Pt ~ Pm ~ P* Pm — 0.27p*; 
thus, we will have a flat universe with a cosmological constant whose density parameter is 
^ 0.73 and, consequently, we can say that we have got a model explaining the same 
observations as the concordance model. In the new model, the dark energy is not vacuum 
energy, but the energy of the A field. Since the p^ value has been fixed, and the relation 
C = p^ is satisfied, Eq. (!2T!) fully defines the function yi = ?/i(C), which can be substituted 
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into Eq. (fT9|l to solve the system formed by this equation and Eq. (fT8|) . In order to find 
a numerical solution of this system, the initial values of y2 and ?/3 are necessary. We have 
taken = 10^ to be well inside the radiation dominated era and after annihilation. As 
it is well known, in this era, the scale factor is proportional to the conformal time; hence, 
one can write a = ar, where a is a constant. Moreover, close to z^, the component Aq will 
be approximately proportional to some power of r, namely, Aq = [3t^ . If these power laws 
are substituted into Eq. ffTOj) . we obtain the relation 5^ — 5 — 6 = 0, whose solution 5^ = +3 
= —2) defines a growing (decaying) mode. We then use the growing mode Aq = j3^T^ 
plus the relation r = (aC)~^ to get the following initial values for the variables j/2, and 

yun = y2^n = yam = -s/J+a-^c^^ . (22) 

Since is constant, its value can be calculated from the initial conditions (!22|) and Eq. (ITSll . 
The resulting formula is 

PA - (23) 

The initial conditions necessary to solve Eqs. (fT8l) and ( |T9i) : namely, quantities i/2m and ysj^ 
are given by Eqs. ( l22l) . but we need the explicit values of the constants a and /?^. In order 
to get these constants, we proceed as follows: since yxm has been already calculated, the first 
of Eqs. (!22|) allows us to obtain a = —yi^^Cir^ then, from the known values of a and p^, 
plus Eq. fl23l) . we get (3^ = ±a^(2p^/25?7)^/^; hence, parameter r] is arbitrary, but it only 
can take on positive values. After fixing 77, there are two /J"*" values with the same absolute 
value and opposite signs. Each of these values generates a set of initial conditions. Indeed, 
in Eq. ( l23l) . it is easily seen that the value of (estimated to explain the observations, see 
above) fixes the product {f3~^)^T], but the rj value and the sign of [3^ remain arbitrary. 

It is evident that, if functions y2(C) and 2/3(C) satisfy Eqs. ffTSj) and f|T9|) . and D is an 
arbitrary constant, functions Dy2{C) and Dy^{Q also satisfy these equations. Moreover, 
according to Eqs. ( !22|) . quantities ?/2m and are proportional to which means that, if 
we find the numerical solution of Eqs. ( fTSil and (fT9l) for /3+ = a^(2p^/25)^''^ (this positive 
Z?"*" value corresponds to 77 = 1), the products of the resulting functions ?/2(C) and yz{C) by 
D = ±l/?7^/^ give other solution. We can thus obtain the solution corresponding to any 77 
value and sign. 

The Aq function corresponding to = a^(2p^/25)^/^ is given in Fig. [1] (in terms of 
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variable (). Equation (fT6|) can be seen as a constraint which must be satisfied by functions 
yi, y2 and y^. We have verified that the functions we have obtained in our numerical 
integration satisfy this equation with a relative error (Si — S2)/S2 smaller than 2 x 10^^^ 
for any z, where Si (H2) is the value of the left- (right-) hand side of Eq. fITB]) . This is a 
satisfactory numerical test for our numerical methods and calculations. Very similar results 
have been obtained with various numerical methods designed to the integration of systems 
of first order differential equations. 



< 




10° 10^ lO'* 10'^ 10*^ 

1 +z 

FIG. 1: Aq component of the four- vector A as a function of 1 + z, for 77 = 1 and positive /?"*". 



V. GENERAL DISCUSSION 

In the concordance model, we have baryonic matter, dark matter, and vacuum energy in 
well-known proportions. The problem with vacuum energy is that theoretical predictions 
— based on standard quantum field theory — lead to very big values of this kind of energy, 
which are not compatible with those required to explain current observation (this is the 
so-called cosmological constant problem). We might imagine some new quantum theory 
leading to a strictly vanishing vacuum energy, but a good enough fitting to the particular 
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value 0.73p* seems to be a fine tuning; however, the vector field is not normalized either 
by theoretical arguments or by non-cosmological observations, which means that, if this field 
exists in nature, cosmological data (p^ value) may be used, as it is done in previous section, 
to succeed in the required normalization. On account of these comments, the new origin for 
the cosmological constant seems to be very good news. 

In this paper we have proved that a certain cosmic field A, with an appropriate La- 
grangian, can play the role of the gravitational constant. This field is coupled to the scalar 
curvature and the Lagrangian have the following form: 

jC = + ^ [R^uA'^A" + V.A^V^A'^] - ^F^uF'^" . (24) 

The constant energy density can take on the right value (0.73p*) whatever the values 
of the rj and e parameters may be. Component Aq evolves (see Fig. [Hand Sec. II VI) . but 

is constant and the relation = —1 is valid at any time. The resulting model is fully 
equivalent to the standard concordance model. 

The new version of the cosmological constant appears in the study of the Universe (cos- 
mology) . In other CcLSGS cLS, e.g., spherically symmetric systems, fully asymmetric structures, 
and so on, the energy-momentum tensor of the field A must be calculated from the formula 
([6]) of Sec. [TTl We have verified that, in the static spherically symmetric case where, in 
adapted coordinates (t, r, 9, 0), the field has the components (y4o(r), Ai(r), 0, 0), the energy- 
momentum tensor does not reduce to the form —p^g^u- Our cosmological constant appears 
in the homogeneous isotropic case; namely, in cosmology. Also in this sense, we are concerned 
with a new type of cosmological constant. 

Perturbations of the field {Aq, 0, 0, 0) are being studied. The vector part of these pertur- 
bations could help to explain the isotropy violations — at large angular scales — detected 
in the WMAP maps of CMB temperatures (see Sec. [T]). In this way, the results of papers 
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12| could be improved. Other nonlinear applications of the theory could lead to bounds, 



relations, or well-defined values of constants rj and e; we are also studying this possibility. 
Finally, generalizations of this theory could be studied, e.g., the field A'^ could be constrained 
to be timelike with the unit norm, and a mass term could be introduced in the Lagrangian 
(see comments in Sec. ilTl). 
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Addendum 

While this paper was in production in Phys. Rev. D, we were advised by A. L. Maroto 
and J. Beltran Jimenez, that similar conclusions were previously obtained by themselves in 
arXiv: 0811.0566vl [astro-ph]. Since the first version of our paper was submitted to Physical 
Review on Nov 30 2008 (26 days after the paper 0811.0566vl was included in the arXiv), 
the publication in PRD has been cancelled. In 0811.0566vl, the cosmological constant was 
interpreted in the framework of a modified theory of the electromagnetic field (Einstein- 
Maxwell generalization), whereas our paper was initially based on an unconstrained vector- 
tensor theory which was not interpreted as a theory of the electromagnetic field. The reader 
can verify that both theories seem to be different; however, we have recently verified -after 
Beltran&Maroto complaints- that they are based on equivalent Lagrangians whose difference 
is a total divergence. Indeed, both papers are based on the general vector-tensor theory 
proposed by C. M. Will several decades ago. We would have never tried any publication 
based on an explicit modification of the Einstein-Maxwell equations. 
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